Disordered systems and the metal-insulator Transition: A super universality class. 



OO 
On 



(N 

c/5 



i 

c 

o 
o 



(N 
> 

o 
o 



c3 



a 

o 

o 



X 
S3 



E. Hofstetter 

Blackett Laboratory, Imperial College, London SW7 2BZ, United Kingdom 

(To be published in Phys. Rev. B) 

The critical behaviour of three-dimensional disordered systems is investigated by analysing the spectral 
fluctuations of the energy spectrum. Our results suggest that the initial symmetries (orthogonal, unitary 
and symplectic) are broken by the disorder at the critical point. The critical behaviour, determinedby 
the symmetry at the critical point, should therefore be independent of the previous invariances and be 
described by a "super" universality class. This result is strongly supported by the fact that we obtain 
the same critical exponent v ~ 1.35 in the three cases: orthogonal, unitary and symplectic. 
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The transport properties in three-dimensional disor- 
dered systems are a very rich topic that already has at- 
tracted considerable research In particular the sys- 
tem exhibits a metal-insulator transition (MIT) as a func- 
tion of the disorder. Although it seems clear now that 
the MIT is a second-order phase transition many points 
still remain unclear. One of them is the influence of the 
symmetry on the MIT. It is generally assumed Q that 
the critical behaviour at the MIT can be cast into three 
different universality classes according to some general 
symmetries of the system: orthogonal (with time reversal 
symmetry, O(N)), unitary (without time reversal sym- 
metry, e.g. with a magnetic field, U(N)) and symplectic 
(with spin orbit-coupling, Sp(N)). One then expects dif- 
ferent critical exponents related to the MIT for the three 
different universality classes. Surprisingly, in spite of the 
change of universality class, the same value of the critical 
exponent has been found, numerically, with and without 
a magnetic field [§,1). The MIT being a transition be- 
tween a chaotic (metallic regime) and a non-chaotic sys- 
tem (insulating regime) a convenient way to study 
this problem is to resort to random matrix theory (RMT) 
and energy level-statistics (ELS). In RMT the statistics 
of the energy spectrum are generally described by three 
different ensembles, Gaussian orthogonal (GOE), unitary 
(GUE) and symplectic (GSE) ensembles depending upon 
the symmetries mentioned above. It has been shown 
HQ , in the case without magnetic field, that besides the 
two expected statistics, namely the GOE for the metallic 
regime and the Poisson ensemble (PE) for the insulat- 
ing regime, a third statistic, called the critical ensemble 
(CE), occurs only exactly at the critical point. The prop- 
erties of the CE will reflect, in particular, the shape of 
the critical states, which are multifractals ||. It was re- 
cently proposed P,[L0j that a natural way to understand 
[^|,D would be that the universality class of the system is 
determined by the symmetry of the system at the criti- 
cal point with the previous results being interpreted 



as the sign of the independence of the universality class 
to the breaking of time-reversal symmetry at the criti- 
cal point. We think that when the system reaches the 
critical point, for increasing disorder, the O(N) or U(7V) 
symmetry is spontaneously broken. The reason for this 
break is obvious when one considers the localised regime. 
If the matrix would be O(N) or U(N) invariant it would 
then be possible to construct delocalised states by lin- 
ear combinations of the localised ones which is clearly 
not allowed. In RMT this phenomenon has led to the 
concept of a preferential basis. The distribution of the 
eigenvalues and the eigenvectors, which are independent 
for the GOE or GUE, become correlated when the disor- 
der is increased and eventually break the O(N) or U(7V) 
invariance [p"l"l| . Similar results concerning the break of 
the O(N) invariance have been reported when studying 
the two-level correlation function Considering the 

spacing distribution P(s) of the energy levels, previous 
results jl(]] suggest that the term inducing the eigenvalue- 
eigenvector correlation and that eventually breaks the 
symmetry appears, at the critical point, only in the large 
s behaviour of P(s) and not in its small s behaviour. 
The symmetry of the system is therefore reflected only 
in its large s behaviour. However, in a recent paper [|13| , 
numerical evidence has been presented showing a small 
difference between the scaling properties of orthogonal 
and unitary systems using more accurate raw data. But 
it has been shown |14j] that an increase of the accuracy of 
the raw data with the transfer matrix method generates 
a systematic shift of the derived quantities and seems to 
give worse results, casting doubts about the real accuracy 
of them. Moreover only one distribution (box), for the 
site energy, was considered in [|l3j whereas three different 
(box, Gaussian and binary) were studied in (TJ]. 

An interesting way to check this problem of symmetry 
is to study and to compare what happens at the critical 
point for the symplectic case or, in other words, to con- 
sider a system with spin-orbit coupling. It is generally as- 
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sumed that the critical exponent for the symplectic case 
should be different (v ~ 1) Jll| from the orthogonal or 
unitary case (v— 1-35) although some recent numerical 
calculations [ |16| have surprisingly given a critical expo- 
nent higher than expected in the symplectic case. In this 
paper we would like to show that the picture with the 
three universality classes to describe the MIT is proba- 
bly inadequate and to propose a new way to understand 
the MIT, shedding new light on the problem. 

In order to investigate the MIT with spin-orbit cou- 
pling we consider the following spin-i tight-binding 
Hamiltonian |17fl: 



ff = 5>.-b>>0>!+ ]T V°f'\j,a)(j',a'\, (1) 



with the hopping matrix elements represented by 



v a f, =t°l 

3,3' 
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where the sites j are distributed regularly in three- 
dimensional (3D) space, e.g. on a simple cubic lattice. 
Only interactions with the nearest neighbours are con- 
sidered. The site energy is described by a stochastic 
variable. In the present investigation we use a box dis- 
tribution with variance = \JW 2 jY2. W represents the 
disorder and is the critical parameter. 1 and er 7 are 
the identity and the Pauli matrices. We choose the cou- 
pling constant (J, = 1, t = 1 and t 1 as independent ran- 
dom variables distributed uniformally between [—5,5] 
for each pair of lattice sites j and j'. The Hermiticity 
and the time reversal symmetry impose Vp^ = {Vy j a )* 

and Vj~p~ a = aa '(^j T f )* respectively, with a, a' = ±1. 

Based on this Hamiltonian, the MIT in the presence of 
spin-orbit coupling is studied by the ELS method, i.e. via 
the fluctuations of the energy spectrum . Starting from 
Eq. (1) the energy spectrum was computed by means of 
the Lanczos algorithm for systems of size M x M x M 
with M = 13, 15, 17, 19 and 21 and disorder W rang- 
ing from 8 to 30. The number of different realizations of 
the random site energies £j j(T was chosen so that about 
10 5 eigenvalues were obtained for every pair of parame- 
ters (M, W) for which only half of the spectrum around 
the band center is considered so that the results do not 
deteriorate due to the strongly localised states near the 
band edges. We checked that we obtain the same re- 
sults using bands of the energy spectrum ranging from 
10 to 50% (half spectrum). After unfolding the spectrum 
obtained, the fluctuations can be appropriately charac- 
terised by means of the spacing distribution P{s) and the 
Dyson-Mehta statistics A 3 (L). P(s) measures the level 
repulsion, it is normalised, as is its first moment, because 
the spectrum is unfolded. A 3 (L) measures the spectral 
rigidity. 

Before studying what happens at the critical point we 
checked, for P(s) and A 3 (L), that one finds, as expected, 



the GSE and the PE regimes for small and large disorder 
respectively. This can be seen in Fig. || for P(s) and Fig. 
I for A 3 (L). 

The next step is to find where the MIT takes place. 
For this, one uses that the quantities we are consid- 
ering here, P(s) and A 3 (L), are scale invariant at the 
critical point. This is due to the fact that the MIT 
is a second order transition and that finite-size scal- 
ing laws apply close to the transition. So we calculate 
a r (M, W) = ± J 3 ° A 3 (L)dL as a function of M and W, 
the critical disorder, W c , being given by W for which 
a r (M, W) is independent of M |Q. a r (M, W) has been 
calculated with a step for W of 0.25 around the critical 
disorder and between 1 and 2 otherwise. But instead of 
using directly these data (raw data) we used the analyt- 
ical properties of a r (M, W) for a finite system to fit a 
third order polynomial between W = 16 and W = 26 
for each M . The results (a(M, W)) are shown in Fig. [|. 
The raw data (a r (M, W)) that where used for the fit are 
given in the inset of Fig. |l|. This method is very inter- 
esting because, as we will see later, less computational 
effort is required and more accurate results are obtained. 
We can already see that the critical disorder W c is accu- 
rately determined. One finds W c = 21.75 ± 0.10. This 
value is higher than in the orthogonal, W c ~ 16.5 [Q, 
and unitary case where it saturates, as a function of the 
strength of the magnetic field, around W c ~ 18.75 [fuj . 
This is due to the weak localization, suppression of the 
weak localization, and antiweak localization phenomena, 
which take place in the orthogonal, unitary, and sym- 
plectic cases respectively jl9| . It is then more difficult to 
localise the states, which means an upwards shift of W c . 
With W c it is now possible to study P(s) at the critical 
point. The first feature we would like to consider is the 
small s behaviour of P(s). In fact for numerical reasons 
what we calculate is the cumulative level-spacing distri- 
bution I(s) = P(s')ds' which allows a better study of 
the small s behaviour of P(s). The results are plotted 
in Fig. |^. We have P(s) oc as 4 with spin-orbit coupling, 
P{s) tx bs 2 with magnetic field and P(s) oc cs without. 
The behaviour is the same as in the metallic regime ex- 
cept for the prefactors a, 6, and c, which are now higher, 
showing a decrease of the level repulsion due to the mul- 
tifractal nature of the critical states The different 
small s behaviours obtained has been interpreted []l6 20 
as a sign that the disorder does not modify the symmetry 
of the system and that we have, indeed, three different 
universality classes at the critical point. We claim here 
that this is not necessarily the case. Considering a sys- 
tem with N sites one sees [pi] that the s@ factor, with 
f3 = 1, 2, 4, responsible for the small s behaviour of P(s) 
is, in fact, a geometrical factor related only to the size of 
the Hamiltonian, N x N, 2N x 2N and 4N x 47V for or- 
thogonal, unitary, and symplectic, respectively. In RMT, 
which corresponds to the metallic regime, this geometric 
factor certainly reflects the symmetry of the system but 
we do not believe it is necessarily the case for a system 
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with large disorder. The distribution of the eigenval- 
ues and the eigenvectors, which are independent for the 
GOE, GUE, or GSE, become correlated when the disor- 
der is increased and eventually break the 0(N), U(A), or 
Sp(iV) invariance. We think that the term inducing the 
eigenvalue-eigenvector correlation and which breaks the 
symmetry at the critical point appears first in the large s 
behaviour of P(s). For larger disorder, in the thermody- 
namic limit, this term will modify the small s behaviour 
of P(s) too, by cancelling the geometrical factor giving 
rise eventually to the PE. So following our argumenta- 
tion what one should compare is the large s behaviour 
of P(s) for the three different ensembles at the critical 
point. The results are reported in Fig. |[ As already 
shown |to| we see no difference between the orthogonal 
and the unitary case. The surprising fact is that we ob- 
tain the same results for the symplectic case too (Fig. ||). 
According to what we wrote above this would mean that 
there is only one "super" universality class at the critical 
point in contradiction to the fact that one expects a dif- 
ferent critical exponent for the symplectic case. A way to 
check that is to calculate v for the symplectic case. Start- 
ing with a(M, W) defined above, it was shown |L8| that 
a can be expressed as a(M,W) = fiM/^^W)) with 
£,oo{W) = Z\W-W C \-", the correlation length. a(M,W), 
being analytical for a finite system, can be written around 
the critical point as 

a(M, W) ~ a(M, W c ) + C\W - W c \M ljv . (3) 

To perform the scaling procedure with the a r (M, W) 
data, the range of a r (M, W) values for various M at any 
given disorder W\ must overlap the range of a r (M, W) 
values for various M for at least one different disorder 
Wi. That would require us to compute a r (M, W) for 
a very large number of W values or for larger M. Us- 
ing the method with the fit as described above this is 
not necessary because from the fit we derived all the val- 
ues we need for the scaling procedure. Using the values 
for a(M, W) given in Fig. [j] we show in Fig. [| that 
a(M, W) can be, indeed, expressed by a scaling func- 
tion /(M/£oo(W)). In the inset of Fig. |^ the correlation 
length £00 (W) is reported as a function of the disorder. 
We clearly see the divergence at W c ~ 21.75. The quality 
of the curves shows that using a fit to the raw data cer- 
tainly improves the results. This is also reflected when 
using Eq. (3) the value of the critical exponent v is de- 
rived. The tricky point is that the formula (3) is valid 
only in the vicinity of the critical point but, on the other 
hand, close to the transition the numerical inaccuracies 
are largest. The task is therefore to find an adequate 
range for \W — W c \ that satisfies these two constraints. 
With the raw data the accuracy of the results as well 
as the value of v is very sensitive to any change in the 
range of \W — W c \. With the value of a(M,W) from 
Fig. [l] this is no longer the case. The quantities calcu- 
lated are much more stable on a wider range of \W — W c \ 
indicating a better accuracy of our results. We obtain 



v = 1.36 ± 0.10 for the critical exponent which in very 
good agreement with v = 1.35 ± 0.15 for the orthogonal 
HI! and v = 1.35±0.20 for the unitary case §. These 
results seem strongly to support that the behaviour of the 
MIT can be cast into one "super" universality class and 
not into three as previously claimed. It has to be noted 
that recently, using a generalized version of the Ando 
model |23[| , a value of v = 1.30 ± 0.20 has been found 
for the symplectic case jl6) in good agreement with our 
results using the Evangelou-Ziman model [0. It is inter- 
esting to note that the results at the critical point seems 
to be independent of the model as already noted in two 
dimensions p4| . 

Concerning A 3 (L) the results, Fig. ||, are more difficult 
to interpret. We know that the shape of Aa(i) contains 
a term linear in L as well as a non-linear term L w , with 
u < 1 ||. As for P(s), it seems we have two different 
behaviours. For the nonlinear term, one obtains differ- 
ent curves for the O(A), U(A) and Sp(A) cases. This 
behaviour was already observed by Batsch et al. |2(J in 
the case with a magnetic field. If now we consider the 
linear term, all the curves seem to fall onto one curve. 
It is interesting to note that the two-point correlation 
function R(s), from which A.$(L) can be obtained, has 
recently been calculated p2| . It was shown that R(s) is 
composed of two terms. One of them is related to the 
break of the O(A) invariance and gives rise to the linear 
term in As(L). So as for P(s) it seems that the break 
of the invariance and therefore the change of symmetry 
is only reflected in one part of As(L), namely the linear 
term. Numerically, even with an accuracy of 0.1 on the 
critical disorder, it is quite difficult to see what happens 
for large L and more work clearly needs to be done to 
increase the accuracy of the results. 

In conclusion we think that the picture, with the three 
universality classes, O(A), U(A) and Sp(A), to describe 
the MIT is probably inadequate. This picture comes from 
a field-theoretical approach using a standard a model 
|l9| , p5| . But we think our results are not necessarily 
in contradiction with the field-theoretical approach. In- 
deed, although these analytical results give indications as 
to the existence of a MIT as well as information about the 
weak localized regime, they say nothing about the criti- 
cal regime (critical point, critical exponent). The prob- 
lem comes from the 2 + e expansion used to solve the a 
model. This perturbative analysis is often unreliable for 
describing the critical regime [E6j even near D = 2. It is 
now well known that the 2 + e expansion gives incorrect 
results for D = 3 when applied to MIT (Ref. |l9|) and 
therefore should be used with great care. Moreover, it 
is far from obvious that this method is adapted for the 
case with magnetic field or spin-orbit coupling where the 
lower critical dimension is < 2. In particular, the multi- 
fractal character of the critical states || , which we think 
plays a crucial role in the description of the critical be- 
haviour, is up to now completely beyond the scope of the 
2 + e expansion. In contrast, interesting progresses have 
recently been made applying supersymmetry techniques 
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to the (non linear)er model to go beyond the perturba- 
tive analysis [^7j although the critical properties of the 
MIT still remain unresolved in the frame of the field the- 
oretical approach J25|,^7]]. Our results suggest that the 
large disorder breaks the O(N), \J(N), or Sp(N) invari- 
ance at the critical point. But this break of symmetry is 
only reflected in some parts of P(s) or As(L). Compar- 
ing these parts one finds that the critical behaviour of the 
MIT is no longer described by three different universality 
classes but by one "super" universality class. This result 
is supported by the fact that one obtains the same criti- 
cal exponent v ~ 1.35 with and without magnetic field as 
well as with spin-orbit coupling. Moreover, Ohtsuki and 
Kawarabayashi recently showed [^8| that the anomalous 
diffusion and the the fractal dimension D(2) is the same 
for O(iV), U(iV), or Sp(N) at the MIT in agreement with 
our results. 

Finally it is interesting to note that our results are in 
good agreement with, at least, some experiments, for the 
value of the critical exponent |^,^0| as well as the the 
absence of influence, at the MIT, of the magnetic field 
pO| , pT| and the spin-orbit coupling p9fl . But clearly fur- 
ther experiments need to be done to check these points 
carefully, in particular to understand the different results 
obtained for uncompensated and compensated semicon- 
ductors. 
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FIG. 1. a as a function of M and W. The critical disor- 
der is given by the point for which a is independent of M. 
W c = 21.75 ± 0.10. The raw data (a r ) used for the fit to 
derive a are given in the inset 

FIG. 2. ln-ln plot of the cumulative level-spacing distribu- 
tion I(s) for small s. We have P(s) oc as 4 with spin-orbit 

coupling (— ■ — ), P(s) tx bs 2 with magnetic field ( ) and 

P(s) ex cs without (••■)• The behaviour is the same as in 
the metallic regime ( — ) except for the prefactors a, b, and c, 
which are now higher. 

FIG. 3. Large s behaviour of P(s). * and + are the curves 
at W — 3 and W = 100, respectively, showing the metallic 
and the insulating regimes for P(s). At the critical point 
the curves are independent of the presence or absence of a 
magnetic field as well as of the spin-orbit coupling. In the 
inset there is a In plot of the tail of P(s). 

FIG. 4. In this figure we show that a(M, W) can be ex- 
pressed by a scaling function f(M/^ ao (W)). In the inset the 
correlation length £oo (W) is reported as a function of the dis- 
order. We clearly see the divergence at W c ^ 21.75. 

FIG. 5. Dyson- Mehta statistics A 3 (L). • is A 3 (L) at 
W = 3 and W = 100. At the critical point one obtains three 
different sets of data which seem to merge into one linear set 
with increasing L. 
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